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ABSTRACT 

We have observed Markarian 421 in January and March 2001 with the 
CANGAROO-II imaging Cherenkov telescope during an extraordinarily high 
state at TeV energies. From 14 hours observations at very large zenith angles, 
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~70°, a signal of 298 ±52 gamma-ray-like events (5.7 a) was detected at £■ > 
10 TeV, where a higher sensitivity is achieved than those of usual observations 
near the zenith, owing to a greatly increased collecting area. Under the assump- 
tion of an intrinsic power-law spectrum, we derived a differential energy spec- 
trum dN/dE = (3.3 ± 0.9stat. ± 0.3,^,*,) x IQ-^^ (^^/iq rY^y^-i^-olU ^i-^i-^o-^syst.) 
ph./cm^/sec/TeV, which is steeper than those previously measured around 
1 TeV, and supports the evidence for a cutoff in the spectrum of Markarian 421. 
However, the 4(7 excess at energies greater than 20 TeV in our data favors a 
cutoff energy of ~8 TeV, at the upper end of the range previously reported from 
measurements at TeV energies. 

Subject headings: BL Lacertae objects: individual (Markarian 421) - gamma 
rays: observations 

1. Introduction 

Markarian 421 (Mrk 421, J1104+3812) is a nearby BL Lacertae object {z = 0.031) and 

was the first extragalactic TeV gamma-ray source discovered (Punch et al. 1992). The TeV 
gamma-ray flux is variable, with flaring behavior observed on time-scales of less than an hour 
(Gaidos et al. 1996). Extensive measurements have been performed by several experimental 
groups based on the imaging Chcrenkov technique (Aharonian ct al. 1999a; Krennrich et 
al. 1999; Piron et al. 2001). Multi- wavelength observations support the Synchrotron-Self- 
Compton (SSC) mechanism for the production of TeV gamma-rays from this source (see, 
e.g., Takahashi et al. 2000; Krawczynski et al. 2001). 

TeV gamma-rays from extra-galactic sources suffer absorption due to photon-photon 
interactions with the inter-galactic infrared background radiation (Nikishov 1962; Gould & 
Schreder 1967; Stecker, de Jager & Salamon 1992). According to recent measurements of the 
infrared background (see, e.g., Hauser & Dwek 2001, and references therein) and predictions 
of the optical depth for TeV gamma-rays (Primack et al. 1999; de Jager & Stecker 2001; 
Totani & Takeuchi 2002), gamma-rays at energies above 10 TeV from Mrk 421 are expected 
to be suppressed, since they interact with mid- to far-infrared photons of ~100 /im. 

Mrk 421 became active in 2000 and 2001, especially at the beginning of 2001 (Borst, 
Getting & Remillard 2001). During this period, northern hemisphere observers measured 

the energy spectrum with good statistics in the region from several hundred GeV to ~10 TeV 
and reported cutoffs at 3-6 TeV (Krennrich ct al. 2001; Aharonian ct al. 2002). The cutoff 
energy is consistent with, or slightly smaller than, that measured for Mrk 501 during its 
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flaring state in 1997 (Aharonian et al. 1999b, 2001). As Mrk 501 has a similar redshift 
{z — 0.034) to Mrk 421, this suggests the cutoffs may be due to infrared absorption of TeV 
gamma-rays. 

We observed Mrk 421 during the 2001 high state with the CANGAROO-II 10 meter 
telescope, at very large zenith angles of ~70°. Similar observations have been reported by 
the Durham group for Mrk 501 in the high state of 1997 (Chadwick et al. 1999). For these 
observations, an effective collecting area ~10 times larger than that for observations near 
the zenith is obtained, with an accompanying increase in the gamma-ray energy threshold 
to -10 TeV. 



2. Observations and Analysis 

The observations were made with the CANGAROO-II 10 meter telescope (Mori et 
al. 2001; Tanimori 2001), located near Woomera, South Australia, Australia (136°47'E, 
31°06'S). The telescope consists of 114 segmented optical mirrors, each of 80 cm diame- 
ter (Kawachi et al. 2001). The camera contains 552 half-inch photomultiplier tubes, arranged 
at 0°.115 intervals, and covering a field of view of —3°. 

Mrk 421 was observed for ten nights in early 2001; January 24, 26, 27, 30, 31 and 
February 1, and March 1-4 (all dates in UT), when the source was extremely active. From 
the CANGAROO-II telescope site, Mrk 421 culminates at a zenith angle of 69°. 3. Approx- 
imately two hours observations were made per night. OFF source data were taken with 
the right ascension suitably offset. An event trigger was registered when 3 individual pixels 
exceeded a threshold of ~2.5 photoelectrons. After rejecting data affected by clouds and 
those at zenith angles greater than 71°. 5, 14.34 hours ON source data and 16.65 hours OFF 
source data remain. A software trigger was applied in order to reduce the effect of pixels 
randomly triggered by the night sky background. Pixels with pulse-heights of greater than 
~3.3 photoelectrons, pixel trigger times within 40 nanoseconds of the central value for the 
event, and three or more adjacent pixels were required. Finally, four or more pixels surviving 
these cuts were required in each event. 

Large zenith angle observations are well-suited to searching for gamma-ray signals at 
higher energies, as a much larger effective area can be achieved compared to observations 
near the zenith (Sommers & Elbert 1987; Tanimori et al. 1994), though with a higher 
energy threshold. From Monte Carlo simulations (Okumura et al. 2001; Enomoto et al. 
2002a), an effective area of ~ 5 x 10^ cm^ at £■ = 20 TeV was estimated for observations at 
70°, with the area increasing to ~ 10^° cm^ for higher energies. A threshold energy (where 
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the gamma-ray detection rate is maximized) of ~11 TeV was derived for a E~^-° spectrum. 
This is an increase by a factor of ~30 in comparison with observations near zenith. The 
energy threshold changes by ~ ±1 TeV if the spectral index is varied by ±0.5. 

The selection of the gamma-ray events is based on the parameterization of the elon- 
gated shape of the Cherenkov light image using the standard parameters: width, length 
(shape), distance (location), asymmetry (direction), and alpha (orientation angle) (Hillas 
1982; Punch 1993; Reynolds et al. 1993). Instead of the conventional parameterization cuts, 
we adopted the Likelihood method (Enomoto et al. 2002a,b), which has a higher efficiency of 
gamma-ray discrimination than the conventional parameterization technique. The likelihood 
method uses a single parameter, Rprob = Prob{'j)/[Prob{'j) + Prob{B.G.)], where Proh{pf) 
and Prob{B.G.) are the probabilities for the event having been initiated by a gamma-ray 
from the source or a background event, respectively. They are the products of individual 
probabilities for width, length, and asymmetry, which are derived from the probability 
density functions, including the energy dependence. These functions were obtained using 
gamma-ray simulations for the signal and the observed OFF source events for the back- 
ground. Rprob ranges from to 1, and the probability of a gamma-ray origin for an event 
increases as Rprob becomes closer to 1, though the gamma-ray acceptance does not signifi- 
cantly decrease in the range of Rprob % 0.5. We adopted a relatively loose cut of Rprob > 0.4 
with an additional requirement of 0°.2 < distance < With these cuts and a further cut 
excluding events with alpha > 20°, 86 % of background events are rejected while 63 % of 
gamma-ray events are expected to be retained. 

The resulting event distribution of alpha is shown in Figure 1 (a) (left panel). A clear 
excess over the background is apparent around the source direction. The excess is broadly 
distributed, up to ~ 30°, due to the deterioration of the pointing resolution, caused by the 
shrinkage of the gamma-ray shower image. This spread in alpha distribution is consistent 
with simulations, as shown in the bottom panel of Figure 1. The OFF source distribution was 
normalized to that of the ON source by the ratio of the number of the events in alpha > 40° 
(0.88), which is consistent with the ratio of observation times (0.86), within statistical errors. 
An excess of 298 ±52 events, with a significance of 5.7 a (calculated using the method 
of Li & Ma (1983)) was obtained in the region of alpha < 20°. For the confirmation 
of the detected signal, the conventional parameterization cuts of 0°.2< distance <l°.l, 
0°.06< length <0°.18, and 0°.03< width <0°.14 were applied to the data, and a signal of 
286 ±55 events was obtained with 5.2 a significance. 

Since the observations were undertaken at large zenith angles, ~ 70°, we carefully ex- 
amined the data and the simulations in more detail: 

1. The shrinkage of the shower image, which is problematic for large angle observations. 
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was studied by a comparison between simulations and data using the background events due 
to cosmic-ray hadrons. Figure 2 shows the imaging parameters length and width, observed 
at large (~70°) and small (~15°) zenith angles, respectively. The hadron simulations were 
made using the CORSIKA code (version 6.004) (Heck et al. 1998), considering the cosmic- 
ray abundance in the TeV region (Mohanty et al. 1998). The resultant distributions of the 
simulations agree with the data for both small and large angles. 

It is also noted that our high resolution imaging camera, which has a pixel spacing size of 
0.115°, helped to separate the smaller images of gamma-ray events from those of background 
events. The expected length and width distributions of the gamma-ray events, simulated 
with the spacing size of 0°.115 and 0°.230, are shown in the bottom panels of Figure 2. With 
the larger pixel size, the reconstructed image size increases and becomes more similar to 
those of hadrons, with an estimated ~50 % decrease in the separation efficiency. Although 
gamma-ray detection is still possible with the larger spacing, the higher resolution imaging 
camera is more advantageous for large zenith angle observations. 

2. The distance distribution of the gamma-ray selected events was compared to those 
from simulations. The location of Cherenkov images due to gamma-ray cascades in the 
field of view has a particular distribution, while those due to hadron showers are uniformly 
distributed. For large zenith angle observations in particular, as the observed distances of 
gamma-ray shower images decrease, there is a substantial difference with the background 
distribution. Figure 3 shows the distance distribution of the gamma-ray selected events, 
which is obtained by subtracting the OFF source distribution from the ON source distribution 
after the likelihood and alpha cuts were applied. The resulting distribution has a clear peak 
around 0°.7 from the source direction, which agrees reasonably well with simulations and 
differs from that of the cosmic ray background, which provides additional confirmation of 
the detection of TeV gamma-rays. 

3. The "standard candle" at TeV energies, the Crab nebula, was observed at relatively 
large zenith angles of ~55° in November and December 2000, and the gamma-ray flux was 
measured for the confirmation of the analysis method and the estimation of the systematic 
error in the energy scale. Using the same analysis technique as that used for Mrk 421, the 
differential energy spectrum was derived over the energy range from 2 TeV to ~20 TeV (Itoh 
et al. 2002), which agrees well with other experiments (Tanimori et al. 1998; Aharonian et 
al. 2000; Krennrich et al. 2001), within a ~15 % error in the energy scale. 

These consistencies provide robust supporting evidence for the detection of > 10 TeV 
gamma-rays from Mrk 421. 
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3. Discussion 

Figure 4 (inserted panel) shows the raw energy spectrum of the observed gamma-ray 
events from Mrk 421. The gamma-ray energy was assigned from the pulse-height sum of the 
individual pixels, using a relation obtained from the simulations. This method is similar to 
that described in Mohanty et al. (1998), and an energy resolution of ~31 % is estimated. 
The excess events are distributed in the energy range 7-45 TeV, however one must take care 
of the spill-over effect from the lower energies due to the finite energy resolution. In order 
to take this effect into account, simulated gamma-ray spectra, with the spectral indexes and 
cutoff energies varied, were compared to the data and the observed spectral parameters were 
determined from the values which minimized the value of x^. With the assumption of a 
power-law spectrum, the differential flux was fitted by 

Jjy / \ -(4.0 Stat. iO.Saj/st.) 

— = (3.3 ± 0.9stat. ± 0.3,,,0 X 10-^3 ( j ph./cmVsec/TeV 

with x^=2.5/2 d.o.f. The cut dependence on Rprob and alpha parameters, and the trigger 
conditions in the simulation, were considered as sources of the systematic uncertainties. The 
systematic errors giving rise to uncertainty in the energy scale such as Cherenkov photon 
scattering in the atmosphere are not included here, but are considered in more detail later. 
The derived spectrum is steeper than those observed at lower TeV energies. The spectral 
shape was tested with a cutoff spectrum of exp(— E/4TeV), as was derived from the 
measurements by the Whipple and HEGRA-CT groups, with the spectral index being the 
hardest one observed during the strong flaring period (Aharonian et al. 2002; Krennrich et al. 
2002). The fitting result did not improve compared to that with the power-law assumption 
(x^=5.0/3 d.o.f.), as an excess of events above 20 TeV is apparent, as shown in Fig 1 (b). 
An excess of 103 ±26 (4.0 a) was observed with alpha < 20°, while 11 events are expected 
for the cutoff spectrum, based on an estimation using the event ratio between 10-20 TeV and 
over 20 TeV. However, if a cutoff energy of 8 TeV is assumed, the consistency with the data 
becomes better (48 events expected for exp(— £'/8TeV)). This cutoff energy is at the 
high end of the range allowed for Mrk 501 (Aharonian et al. 1999b, see also Aharonian et al. 
2001). Since these two AGNs have similar redshifts, the cutoff energies in both spectra are 
expected to be similar, assuming the attenuation is predominantly due to infrared absorption. 
As there is only a 2 cr difference between our observations and this prediction, our result falls 
in the acceptable range of the absorption hypothesis due to the cosmic infrared background. 

Figure 4 (main panel) shows the measured energy flux, assuming the power-law spec- 
trum. Data for the Whipple (Krennrich et al. 2001) and HEGRA-CT groups (Aharonian 
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et al. 2002), observed during a similar period of the flaring state (January-March 2001) 
are also shown. The observation periods were not exactly the same and the source varied 
significantly during this high state, therefore the absolute fluxes are expected to differ at 
some level. The absolute flux level determined from the CANGAROO-II data is within the 
observed range of the flux variation reported by the Whipple group (Krennrich et al. 2002), 
and the spectral slope around 10 TeV is consistent with that of these two groups, supporting 
the roll-over from the flatter spectrum measured at lower energies. 

For large zenith angle observations, a large uncertainty in the energy scale, due to the 
absorption of Chcrcnkov photons in the atmosphere, is inevitable. Only Rayleigh scattering 
was considered in the simulation code to avoid over-estimating the gamma-ray energies. The 
inclusion of Mie scattering and ozone absorption would affect the energy scale by ~30 % and 
~3 %, respectively, based on numerical estimations using the program code of Kneizys et al. 
(1996). We stress that these effects increase the energy scale. The use of the "flat-Earth" 
approximation for the atmosphere in the simulations requires a ~6 % correction which has 
already been taken into account in the discussion above. 

The measurement of spectra at large zenith angles was verifled by observations of the 
Crab nebula up to the zenith angles of ~55°, although calibration using the Crab nebula at 

the same zenith angles as the Mrk 421 observations (~70°) is unfortunately impractical with 
the current instrumental sensitivity. The strong gamma-ray emission of Mrk 421 (~ 3 times 
that of Crab nebula) enabled us to detect the source in only 14 hours. In order to detect 
the Crab nebula at the same signiflcance level, more than 150 hours observations would be 
required. 

In summary, owing to the large effective area and the high resolution performance of 
the Cherenkov imaging camera, E>10 TeV gamma-rays from Mrk 421 were detected at a 
high confidence level at zenith angles of ~70° with 14 hours of observations. The derived 
spectrum in the region of 10-30 TeV is steeper than that around 1 TeV, which supports the 
cutoff spectrum of Mrk 421 measured in the 0.2-10 TeV range by other groups. The excess 
observed above 20 TeV is strongly suggestive of a higher cutoff energy, ~8 TeV, compared 
to the lower energy observations. These observations confirm, with the support of detailed 
simulations, the viability of the large zenith angle technique. Large zenith angle observations 
provide a unique method of measuring the spectrum in the important energy range above 
10 TeV with a relatively short observation time. 
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Fig. 1. — Image orientation angle (alpha) distributions for gamma-ray-like events with 
respect to the direction to Markarian 421. The left figure (a) shows the distributions for all 
energies, and the right figure (b) for those with reconstructed energies above 20 TeV. In the 
upper panel, filled circles with error bars (statistical only) and solid lines are for the ON and 

OFF source data, respectively. The lower panel shows the excess events of the ON source 
above the background (OFF source) level. The solid curves show the expected spread of 
gamma-ray events in the alpha distribution from simulations. 
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Fig. 2. — Distributions of the image-centroid parameters {length and width) observed at the 
large (~70°) and small (~15°) zenith angles. In the upper and middle panels, the observed 
background data (dots with error bars) and hadrons simulations (solid lines) are shown for 
the large and small angles. In the bottom panels, those of the gamma-ray simulations, with 
different camera pixel spacings (0°.115 and 0°.230), are shown for large angle observations. 
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Fig. 3. — Distributions of the image shape parameter distance, after subtracting normahzed 
OFF-source data from ON-sourcc data (circles with error bars), gamma-ray simulation (solid 
line) and hadron background (dotted line). 
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Fig. 4. — The observed gamma-ray fluxes (main panel) and the energy spectrum of gamma- 
ray events (inserted panel). In the inserted panel, data are represented by circles with error 
bars, with a 2 cr upper limit plotted at the highest energy. Best-fit spectra for a power- 
law (£^~^ ° ; dot-dashed hne) and a cut-off {E"^-^ ex-p{—E/4 TeV) ; dotted line) are shown 
(see text for details). The data shown with the filled circles were used for the spectral 
shape fitting. In the main panel, the measured flux under the assumption of a power-law 
spectrum is shown with error bars and the area corresponding to statistical errors of ±1 a. 
Whipple (Krennrich et al. 2001) and HEGRA-CT (Aharonian et al. 2002) spectra measured 
in similar periods are also shown. The fluxes plotted for the HEGRA-CT group have been 
scaled in order to normalize it to the Whipple flux at 1 TeV. 



